%%  Afgeleide van w
function [f]=func2U8(y,const)
DA=const.mu*y(1)*y(2)-const.k_d*y(1)-const.q*y(1);
DP=-const.mu*y(1)*y(2)+const.q*(const.P_in-y(2));
DD=const.k_d*y(1)-const.k_s*y(3)-const.q*y(3);

f=[DA,DP,DD]';
end